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ABSTRACT 

We present deep 1.8 cm (16 GHz) radio continuum imaging of seven young stellar objects in 
the Taurus molecular cloud. These objects have previously been extensively studied in the 
sub-mm to NIR range and their SEDs modelled to provide reliable physical and geometrical 
parametres. We use this new data to constrain the properties of the long-wavelength tail of the 
greybody spectrum, which is expected to be dominated by emission from large dust grains in 
the protostellar disk. We find spectra consistent with the opacity indices expected for such a 
population, with an average opacity index of P — 0.26 ±0.22 indicating grain growth within 
the disks. We use spectra fitted jointly to radio and sub-mm data to separate the contributions 
from thermal dust and radio emission at 1.8 cm and derive disk masses directly from the 
cm-wave dust contribution. We find that disk masses derived from these flux densities under 
assumptions consistent with the literature are systematically higher than those calculated from 
sub-mm data, and meet the criteria for giant planet formation in a number of cases. 
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1 INTRODUCTION 

The growth of dust grains from sub-micron to micron sizes in the 
local environment of protostellar objects is well-probed by obser- 
vational data in the mid-infrared, however the further growth of 
these grains to larger sizes is only probed by their sub-mm and cm- 
wave emission. It is this growth, and the possible coagulation and 
sedimentation that follows it, which is widely thought to lead to 
the formation of proto-planetary objects in the dense disks around 
evolving protostars. Measuring the long wavelength tail of the ther- 
mal emission from dust, which is expected to be dominated by such 

* We request that any reference to this paper cites "AMI Consortium: Scaife 
etal. 2011". 
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large dust grains, often referred to as pebbles, provides an indepen- 
dent measure of the dust composition in the disk through its opacity 
index (/3), as well as the mass of the accumulated matter in the disk 
which is an important factor in models of planet formation (e.g. 
Boss 1998). 

Sub-mm measurements towards circumstellar disks are often 
used to determine their potential for planet formation, as longer 
wavelength information provides a unique perspective on the cooler 
dust in the outer part of the disk. It is this region where protoplanets 
are expected to form. Alternative methods of disk mass estimation, 
such as spectral line measurements of molecular gas, are compli- 
cated by opacity effects (e.g. Beckwith & Sargent 1993) and re- 
quire complex models to account for these effects as well as those 
of depletion (Dutrey 2003; Kamp & Dullemond 2004). Sub-mm 
and radio measurements are useful not only as a probe of the disk 



2 Scaife et al. 



mass itself, but with multi-wavelength data available, they can also 
be used to determine the evolution of the opacity index as a func- 
tion of frequency (Shirley et al. 2011a;b) and place constraints on 
the growth of dust grains in such disks. This can be used to deter- 
mine whether the timescales assumed in current models of planet 
formation are consistent with observational data. However, to date 
the disk masses determined from sub-mm data appear to be too 
low to agree with theoretical models of planet formation (see e.g. 
Andrews & Williams 2005), although there are a number of obser- 
vational caveats associated with these data. 

In practice, multiple frequencies are required to constrain the 
properties of dust emission and draw comparisons with models of 
radiative transfer (Shirley et al. 201 la). In the sub-mm the observa- 
tional properties of the dust derived from its spectral energy distri- 
bution (SED) can have contributions from both the envelope and the 
disk, which need to be separated. In order to probe the properties of 
the disk alone it is necessary to go to longer wavelengths in order to 
ensure that it provides the dominant contribution to the SED. How- 
ever, this transition regime, between the cm- and sub-mm- wave- 
lengths, is often confused by the presence of multiple emission 
processes each contributing to the overall SED. For younger pro- 
tostellar objects, such as Class sources, (partially-) optically thin 
free-free emission is frequently observed to dominate the cm-wave 
spectrum; for Class I and II objects the high magnetic field strength 
on the surface of the central star can lead to gyrosynchrotron emis- 
sion from quasi-relativistic or thermal plasma. 

In this work we present SEDs for eight Class I protostars in the 
Taurus molecular cloud. These objects were selected from the sam- 
ple of Gramajo et al. (2010) who used multi-frequency data from 
mid-infrared (MIR) to sub-mm wavelengths to model the SEDs of 
these objects. We here extend these SEDs to the cm-wave regime 
with the AMI-LA in order to constrain the long-wavelength prop- 
erties of the optically thin emission expected to be dominated by 
their circumstellar disks. 



2 OBSERVATIONS 

AMI comprises two synthesis arrays, one of ten 3.7 m antennas 
(S A) and one of eight 13 m antennas (LA), both sited at the MuUard 
Radio Astronomy Observatory at Lord's Bridge, Cambridge (AMI 
Consortium: Zwart et al. 2008). The telescope observes in the band 
13.5-17.9 GHz with eight 0.75 GHz bandwidth channels. In prac- 
tice, the two lowest frequency channels (1 & 2) are not generally 
used due to a lower response in this frequency range and interfer- 
ence from geostationary satellites. The data in this paper were taken 
with the AMI Large Array (AMI-LA). 

AMI-LA data reduction has been described extensively in pre- 
vious works (see e.g. AMI Consortium: Zwart et al. 2009). For 
this work flux (primary) calibration is performed using short ob- 
servations of 3C48 and 3C147. From other measurements, we find 
the flux calibration is accurate to better than 5 per cent (AMI 
Consortium: Scaife et al. 2008; AMI Consortivun: Hurley-Walker 
et al. 2009). Phase (secondary) calibration is carried out using in- 
terleaved observations of strong point sources, see Table 1. After 
calibration, the phase is generally stable to 5° for channels 4-7, 
and 10° for channels 3 and 8. In the following analysis chaimels 
4 — 8 have been used. We have applied an increased absolute cal- 
ibration error of 10 per cent to flux densities from channel 8 in 
order to reflect the poorer phase stability of this channel relative to 
the others. 

Reduced data were imaged using the AIPS data package, and 



CLEAN deconvolution was performed using the task IMAGR. The 
half power point of the AMI-LA primary beam is fs; 6 arcmin at 
16 GHz, and in what follows we use the convention: Sy °= v", 
where 5v is flux density, V is frequency and a is the spectral in- 
dex. All errors quoted are 1 a. 



3 RESULTS 

Flux densities were recovered from both the combined frequency 
data and the individual frequency channels. The combined fre- 
quency image is made in the same manner as the individual chan- 
nel images, from the weighted combination of multi-channel vis- 
ibility data. Star-forming regions such as the Taurus molecular 
cloud often display complex emission structures in the radio and 
consequently the data were mapped to identify possible confusion 
with neighbouring sources. In the cases where targets are point- 
like and isolated (IRAS 04016-1-2610, 04361-1-2547, 04302-1-2247 & 
04325-1-2402) flux density estimates from the maps were clarified 
by direct fits to the visibility data; these values were found to devi- 
ate by < 5 per cent. 

Six of the seven target sources were detected in the continuvun 
data, the exception being IRAS 040248-1-2612. IRAS 043024-2247 
was detected in the continuum data, but at too low significance to 
recover reliable detections at > 5 cr in the individual channels. The 
remaining five objects were all detected at > 5 C7 in at least four 
channels, see Table 2. 

3.1 Spectral Indices and Dust Disk Masses 

Where possible, power-law spectral indices were fitted to the AMI- 
LA channel data for each object, see Table 2, using the MCMC 
based Maximum Likelihood algorithm METRO (Hobson & Bald- 
win 2004). The resulting spectra were all found to be steeply ris- 
ing with frequency, with a weighted average spectral index a = 
1.30 ±0.63. Such steeply rising spectra are compatible with free- 
free emission, which has a maximum slope of a = 2 in the optically 
thick regime. The presence of purely optically thick free-free emis- 
sion across the AMI-LA band would imply the presence of ultra- or 
hyper-compact Hii, for which a spectral index of a I is typical. 
However, both of these are more usually associated with the later 
stages of more massive star formation. 

Spectral indices with a > 2, such as that measured for 
IRAS04361-I-2546, are consistent with the slope of the SED pro- 
duced by thermal emission from dust grains. In the Rayleigh-Ieans 
region the opacity of the thermal emission from dust grains can be 
well approximated by a power law, KyocyP. This index p is related 
to the spectral index of flux density measurements, here defined as 
5v °= v", as j3 ~ (1 + A) X (a — 2), where A is the ratio of optically 
thick to optically thin emission. This ratio has been previously de- 
termined in the region 350 /xm to 1.3 mm as A ~ 0.2 (Rodmann 
et al. 2006; Lommen et al. 2007); however, we here assume that 
for much longer cm-wavelengths A — > 0, and that the emission is 
entirely optically thin. If the emission at 1.8 cm is entirely due to 
thermal dust then the measured a would imply values of jS 0. 
Although negative values of j8 are unphysical, values approach- 
ing zero are proposed for the long wavelength tail of the greybody 
spectrum from young stellar objects (e.g. Pollack et al. 1994). Such 
values are low compared with the canonical value for the interstel- 
lar medium of jSism ~ 1.8, due to the assumed change in response 
to the grain size distribution within the dust population, with larger 
dust grains producing shallower opacity indices. Indeed for grains 
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Table 1. AMI-LA sample of Class I protostars in Taurus. Column [1] source name; column [2] alternative designations for the source; column [3] source Right 
Ascension; column [4] source Declination; column [5] major axis of the AMI-LA PSF in the combined channel image; column [6] minor axis of the AMI-LA 
PSF in the combined channel image; column [7] noise on the combined channel image; column [8] date of observation; column [9] primary (flux) calibrator 
source; & column [10] secondary (phase) calibrator source. 



Name 


Other 
names 


RA 
(J2000) 


Dec 
(J2000) 


(arcsec) 


^min 

(arcsec) 


^^rms 


Date 
dd-mm-yy 


Flux 
Cal. 


Phase 
Cal. 


IRAS04016+2610 


L1489 IRS 


04 04 42.9 


+26 19 02 


38.9 


24.6 


23 


26-03-11 


3C147 


J0406+2728 


DGTauB 




04 27 02.6 


+26 05 31 


32.9 


21.6 


27 


07-07-11 


3C48 


J0429+2724 


IRAS04248+2612 


HH31 


04 27 57.3 


+26 19 19 


33.7 


22.1 


25 


28-03-11 


3C147 


J0429+2724 


IRAS04302+2247 




04 33 16.5 


+22 53 21 


33.5 


21.5 


32 


04-04-11 


3C147 


J0438+2153 


IRAS04325+2402 


L1535 IRS 


04 35 35.4 


+24 08 19 


33.4 


23.1 


22 


15-04-11 


3C48 


J0438+2153 


IRAS0436 1+2547 


TMR 1 


04 39 13.9 


+25 53 21 


33.6 


22.2 


30 


18-04-11 


3C48 


J0440+2728 


IRAS04368+2557 


L1527 


04 39 53.6 


+26 03 06 


32.2 


22.6 


30 


19-04-11 


3C48 


J0440+2728 



Table 2. Measured and derived properties of the AMI-LA Class I Taurus Sample. Column [1] source name; columns [2-6] measured flux densities at cm- 
wavelengths; column [7] flux density from the AMI-LA combined channel data; column [8] spectral index across the AMI-LA band. 



Frcq. IGHz] 



Name 


8.44" 


14.62 


15.37 


16.12 


16.87 


17.63 


Sl6 


0!AMI 




(MJy) 


(MJy) 


(MJy) 




(Hiy) 


(MJy) 


(MJy) 




IRAS04016+2610 


343 ±33 


528 ±57 


555 ±57 


510±46 


563 ±44 


615 ±65 


546 ±36 


0.65 ±0.27 


DGTauB 




523 ±49 


555 ±45 


595 ±52 


645 ±60 


735 ±78 


647 ±42 


1.79 ±0.27 


IRAS04248+2612 


< 100 












< 125 




1RAS04302+2247 


160±33 












249 ±28 




1RAS04325+2402 


< 100 


148 ±52 


456 ±50 


299 ±44 


327 ±57 


<345 


338 ±28 


1.49±1.10 


1RAS0436 1+2547 


140 ±32 


382 ±49 


435 ± 36 


504 ±39 




690 ±90 


449 ±37 


2.82 ±0.72 


1RAS04368+2557 


8 10 ±30'' 


918±62 


1040 ±62 


1194 ±60 


1014±71 


1151 ±130 


904 ±54 


1.17 ±0.42 



" 8.44 GHz flux densities from Lucas, Blundell & Roche (2005). 
* 8.5 GHz flux density from MeUs et al. (2011). 



with a size distribution An/Aa oc a^P , where the grain size a ^ Omax 
and Oniax J? 3A, with A the wavelength of observation, the measured 
opacity index, J3 (A) ~ — 3)J3ism- For a tjfpical power-law index 
of p = 3.5, grains with jSism ^ 2 will have j3 «i 1 when Cmax > A; 
and by extension jS < 1 when Omax > 3A (Draine 2006). Conse- 
quently, opacity indices of J3 ^ 1 indicate the largest grains, which 
are expected to exist as a consequence of grain growth in protoplan- 
etary disks. Furthermore, power-law indices for grain size distribu- 
tions as low as p = 3 have been determined in protoplanetary disks 
(e.g. Tanaka et al. 2005) reducing the expected observable value 
of j3 even further. Indeed, Andrews & Williams (2005) determined 
an average jS = 0.06 ± 0.02 in the wavelength range 450-850 jUm 
consistent with these predictions, and jS < in the range 850 ^im 
to 1.3 nmi, the unphysical nature of which they attributed to con- 
tamination of the 1.3 mm data by alternative radio emission mech- 
anisms; a possibility made likely by the spread in protostellar class 
across their much larger sample. 

3.2 Spectral Energy Distributions 

For each object in this sample we use data at 12 ^ A ^ 1300 /im 
from Gramajo et al. (2010) and the SED fitting tool of Robitaille 
et al. (2007) to fit the sub-mm SED. This tool fits SEDs based on 
a gridded procedure with the option of 200,000 SED models, 
composed of 20,000 radiative transfer models (Whitney et al. 2003) 
viewed at a range of 10 inclination angles. Each SED has contribu- 
tions from a number of components including the central star itself, 
the envelope and the disk, the contribution from each of which will 
vary according to evolutionary stage as indicated by the data. 



Since these SEDs are restricted to A ^ 1300 /im, we make the 
assumption that the dust SED is dominated at long wavelengths by 
the contribution from large grains predominantly found in the disk. 
We therefore extrapolate the disk component of the best-fitting Ro- 
bitaille et al. (2007) SED as a power law with a spectral index 
denoted a' (to distinguish it from the low frequency power law 
index, a) from A = 1300/im to the AMI passband to investigate 
whether the AMI-LA data is consistent with the tail of the grey- 
body spectrum. Under an assumption of optically thin emission at 
cm-wavelengths, this index is assumed to be related to the opacity 
index as j8' = a' — 2. In summary, data in the centimetre to sub-mm 
wavelength range are consequently fitted using a double power-law 
with low frequency data dominated by emission with index a and 
higher frequency emission dominated by emission with index a' . 
The resulting SEDs are discussed for each object in the following 
section. 



3.2.1 IRAS 04016+2610 

Lucas, Blundell & Roche (2000; hereinafter LBR2000) previously 
detected IRAS 04016+2610 at 3.5 and 2 cm, and placed a 3 cr 
upper limit on the 6 cm flux density from IRAS 04016+2610 of 
'S4.89GHZ < 260 ^Jy. However, they qualified this non-detection as 
an effect of the high spatial resolution of their MERLIN observa- 
tions as demonstrated by a lower spatial resolution detection by 
Rodriguez et al. (1989) with 54.89GHz = SOOiSO^uIy. The mea- 
sured flux densities from the detections at these frequencies im- 
ply a radio spectrum with a| gg = —0.69 ± 0.70, suggesting a non- 
thermal emission mechanism. The 2 cm flux density of 5i5ghz = 
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IRAS 04016+2610: 



DG Tau B: 




ICO 1000 10* 1 10 100 
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(a) (b) 

Figure 1. Top: AMI-LA 1.8 cm maps of (a) IRAS 04016+2610 & (b) DG Tau B. Contours are shown in increments of fj from 5 a, where a for each object is 
listed in Table 1. The position of the target object is indicated by a white cross ('+'), the FWHM of the AMI-LA primary beam by a circle, and the synthesized 
beam is shown as a filled ellipse in the bottom left comer of each map. These maps are uncorrected for primary beam attenuation. Bottom: Spectral energy 
distributions for (a) IRAS 04016+2610 & (b) DG Tau B. Data points at V > 200GHz are taken from Gramajo et al. (2010), data points at lower frequencies 
are compiled from the literature (see text for details) and this work, see Table 2, with data from the AMI-LA indicated in red in the online version. A best 
fitting SED from the grid models of Robitaille et al. (2007) is shown as a sohd line at v > 200 GHz, with the disk component of this model shown as a dashed 
line. Below V = 300 GHz this disk component is extrapolated as a power-law, see text for details, and shown as a dot-dash line. A second power-law spectrum 
is shown fitted to lower frequency radio data at V < 10 GHz with a = — 1.1 and a = —0.1 for IRAS 04016+2610 and DGTau B, respectively. The combined 
spectrum from the two power-law components is shown as a solid line at V < 300 GHz. 



700 ± 100 /xJy from Rodriguez et al. (1989) is broadly consis- 
tent with that measured by the AMI-LA, and that of LBR2000: 
•^iSGHz = 520±40/xJy. The wider cm-wave spectrum incorporat- 
ing all these data, see Fig. 1(a), may be fitted assuming a domi- 
nant contribution to the SED with a < — 1.1 atv< 10 GHz, above 
which frequency the greybody tail from thermal dust emission be- 
comes dominant. We find a best fitting extrapolation to the disk 
SED with a' = 2.4 (j8' = 0.4), see Fig. 1(a). At face value, these 
spectral indices imply an SED where a non-thermal emission mech- 
anism is dominant at low radio frequencies. 

Synchrotron emission would be unusual in the vicinity of pro- 
tostellar objects as it requires both high magnetic field strengths 
and a population of ultra-relativistic particles. A spectral index 



of 0! = — 1.1 is broadly consistent with gyrosynchrotron emission 
from a quasi-relativistic power-law population of electron plasma, 
which is expected to exhibit steeply falling spectral indices (Dulk 
1985). Gyrosynchrotron has been invoked as a potential radio emis- 
sion mechanism for a number of protostellar objects (e.g. Ray et al. 
1996). The position of the peak in the spectrum of such emission is 
a strong function of the magnetic field strength. A spectrum which 
peaks below 10 GHz, such as that of IRAS 04016+2610, is likely to 
arise from a magnetic field with a magnitude ~ 10 G. The magnetic 
field strength on the surface of Class 1 objects is thought to be of 
order 10^ G, as seen towards the Class 1 object WL 17 (Johns-Krull 
et al. 2009); a field strength two orders of magnitude smaller is 
expected from the inner radius of the circumstellar disk (Donati & 
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IRAS 04248+2612: IRAS 04302+2247: 




Frequency [GHz] Frequency [GHz] 

(a) (b) 

Figure 2. Top: AMI-LA 1.8 cm maps of (a) IRAS 04248+2612 & (b) IRAS 04302+2247. Contours are shown in increments of a from 5 fJ, where a for 
each object is listed in Table 1. The position of the target object is indicated by a white cross ('+'), the FWHM of the AMI-LA primary beam by a circle, and 
the synthesized beam is shown as a filled ellipse in the bottom left comer of each map. These maps are uncorrected for primary beam attenuation. Bottom: 
Spectral energy distributions for (a) IRAS 04248+2612 & (b) IRAS 04302+2247. Data points at V > 200 GHz are taken from Gramajo et al. (2010), data 
points at lower frequencies are compiled from the literature (see text for details) and this worlc, see Table 2, with data from the AMI-LA indicated in red in 
the online version. A best fitting SED from the grid models of Robitaille et al. (2007) is shown as a solid line at v > 200 GHz, with the disk component of this 
model shown as a dashed Une. Below V = 300 GHz this disk component is extrapolated as a power-law, see text for details, and shown as a dot-dash line. A 
second power-law spectrum is shown fitted to lower frequency radio data at V < 10 GHz with a = —0.1 for IRAS 04302+2247. The combined spectrum from 
the two power-law components is shown as a solid line at V < 300 GHz. 



Landstreet 2009) and this may be the source of the gyrosynchrotron 
emission suggested here. 



3.2.2 DG Tail B 

A variable free-free component is suggested by the SED of 
DG Tau B, see Fig. 1(b). The magnitude of the emission measured 
with the AMI-LA is only ~70 per cent of that measured previously 
by Rodriguez et al. (1996), who found 5i5ghz = 0.8 ±0.1 mJy for 
this source, see Fig. 1(b), although the data are consistent within 
3 CT. Variability is also evident in the spectrum of the nearby source 
DG Tau A, which lies immediately to the north-east of DG Tau B. 



The spectrum for this object is presented in Fig. 3, where the spread 
in measurements at 5 GHz illustrates the degree of variability. 

The sub-mm to millimetre SED of DG Tau B is poorly con- 
strained due to its proximity to DG Tau A, however we find that 
the spectrum at frequencies v > 20 GHz can be well described by 
a power-law with a — 2.0(/3 ~ 0.0). This would imply a disk com- 
ponent with emission approximately one order of magnitude higher 
than that of the best-fitting grid SED found by Gramajo et al. (2010) 
and that the disk contribution is significant in DG Tau B. We illus- 
trate this scaled disk model as a dashed line in Fig. 1(b). 

Unlike the Class I sources targeted by this sample DG Tau A 
can be seen to possess a more shallowly rising spectrum in the cm 
to mm-wave regime. This spectrum is nominally well fitted by an 
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1 10 100 1000 10* 10^ 10^ 



Frequency [GHz] 

Figure 3. Spectral energy distribution of DG Tau A. Millimetre and sub- 
millimetre data are taken from the literature, see text for details, archival 
centimetre radio data from Wendker et al. (1995) are shown, with data 
points from the AMI-LA (this work) indicated in red in the online version. 
A best fitting model from the Robitaille et al. (2007) grid is shown as a solid 
line with the disk component indicated as a dashed line. Below 300 GHz this 
disk component extrapolated as a power-law with a' = 2.5(j3' = 0.5) and 
a free-free component with a = 0.6 are shown as dot-dash Hnes. The com- 
bined spectrum from the two power-law components is shown as a solid 
line. 



index of a = 0.6, see Fig. 3, the canonical value for partially opti- 
cally thick free-free emission from a spherical stellar wind (Pana- 
gia & Felli 1975; Reynolds 1986). DG Tau A is a Class II source 
with a well known wide-tailed molecular outflow, thought to host a 
large disk wind (Agra-Agoade et al. 2009), which may be provid- 
ing the source of this emission. Radio flux densities for this source 
are taken from the catalogue of Wendker (1995). We do not include 
DGTau A in our wider analysis as it is not a member of the original 
sample. 



3.2.3 IRAS 04248+2612 

IRAS 04248-1-2612 is a non-detection at 1.8 cm, which constrains 
an extrapolation of the greybody spectrum to having a spectral in- 
dex a' > 2.7 (j3' > 0.7), see Fig. 2(a). 



3.2.4 IRAS 04302+2247 

IRAS 04302-1-2247 is detected at > 5 CT in the AMI-LA continuum 
data only and it is therefore not possible to measure a spectral in- 
dex at cm-wavelengths directly from the AMI-LA data. Previously 
it was also detected at 8.44 GHz by LBR2000, who did not detect 
it at 2 or 6 cm. Assuming that the detection at 8.44 GHz represents 
the contribution of an optically thin free-free emission component 
we fit the spectrum shown in Fig. 2(b). Assuming a canonical opti- 
cally thin spectral index of a = —0. 1 for the free-free component, 
this model is compatible with the AMI-LA continuum data point 
assuming an extrapolation of the best fitting disk component to the 
SED with a spectral index a' = 2.3 (jS' = 0.3). However we note 
that the low frequency spectral parametres for this object remain 
largely unconstrained by these data. 



3.2.5 IRAS 04325+2402 

The AMI-LA image of IRAS 04325-1-2402 is dynamic range limited 
by the bright radio source at the half power point of the primary 
beam in the north-east of this field, see Fig. 4(a). We tentatively 
associate this source with the emission line star Coku-Tau 3. As a 
result the data, and hence image quality per channel, for this source 
are comparatively poorer than remainder of the sample and there 
is increased scatter in the cm-wave spectrum. This is reflected in 
a much larger uncertainty on the cm-wave spectral index. We find 
that the cm-wave flux density can be accounted for entirely by an 
extrapolation of the best-fitting disk component from the Robitaille 
et al. (2007) model with a' = 2.1 (jS' = 0.1). This conclusion is 
also supported by a non-detection of this object at 8.44 GHz by 
LRB2000. 

3.2.6 IRAS 04361+2547 

IRAS 04361-1-2547, see Fig. 2, is found to have a steeply rising 
spectrum across the AMI band with Wami = 2.82 ±0.72. Such a 
spectrum is too steep to be attributed to free-free emission, however 
it is similar to that found for IRAS 16293-2422B (Rodriguez et al. 
2005) 

IRAS 04361-1-2547 was previously detected in the radio at 
4.89 GHz (Terebey et al. 1992) and at 8.44 GHz by LBR2000 
who also made a tentative detection at 2 cm. With a measured 
flux density of 54.89ghz = 230±36MJy (Terebey et al. 1992), a 
steeply falling spectral index is found below 10 GHz with ai gg = 
— 0.92±0.47. Such a spectrum implies a non-thermal emission pro- 
cess such as synchrotron emission, however this would be unusual 
for a protostellar object. Indeed, when considered as part of the 
wider cm-wave SED, for this emission to be compatible with that 
seen at 16 GHz by the AMI-LA it must possess a spectrum with 
a < —3.5, see Fig. 4(b). As described previously, such a spectral 
index would be typical of gryosynchrotron emission from a power- 
law population of electrons; it is also not incompatible with that of 
thermal gyrosynchrotron emission, which can have an even more 
steeply inverted spectrum with ajH = ^8. 

The inference of a non-thermal emission mechanism from the 
SED does not take account of any possible variable free-free emis- 
sion from these objects. The observations at 4.89 and 8.44 GHz 
were taken 7 years apart, and over a decade prior to those of the 
AMI-LA. Variable radio emission has been observed from a va- 
riety of protostellar objects (e.g. Shirley et al. 2006) and it is not 
certain that spectral behaviour determined from non-simultaneous 
measurements in this region of the spectrum can be relied upon. We 
note that the tentative detection at 2 cm of IRAS 04361-1-2547 by 
LBR2000 is broadly consistent with the AMI-LA flux densities (al- 
beit at low significance, < 4.5 CT), and was made at a similar epoch 
to that at 8.44 GHz which can also be accounted for by an extrapo- 
lation of the AMI-LA power-law. This may imply that the spectral 
behaviour of IRAS 04361-1-2547 was similar at the LBR2000 epoch 
to that of the AMI-LA, but different at the epoch of the 4.89 GHz 
observations. Data at frequencies > 5 GHz can be accounted for 
by an extrapolation of the greybody tail using a power-law with 
a' = 2.0(/3' = 0.0), and this is illustrated in Fig. 4(b). 

3.2.7 IRAS 04368+2557 

IRAS 04368-1-2557 is defined by Gramajo et al. (2010) as a bor- 
derline Class O/I object and this is evident from the enhanced en- 
velope contribution to its SED compared with the other sources in 



this sample, see Fig. 5. This object is also part of the EVLA sur- 
vey of edge-on protoplanetary disks (Melis et al. 201 1), which pro- 
vides flux density measurements for this source at 5, 8.5, 22.5 and 
43.3 GHz. These data are combined with the new AMI-LA mea- 
surements presented here and are shown in Fig. 5. Melis et al. 
(2011) fitted data at v < 8.5 GHz and > 8.5 GHz separately and 
constrained the low frequency radio emission to have a spectral in- 
dex of a|-^ = 0.33 ± 0. 17 and the high frequency emission to have 
an index of 0!v>8.5 = 2.87 ±0.17. In this work we fit data from 
5 GHz to 350 GHz jointly using the sum of two power-laws. There 
is a degeneracy between the fitted indices and, using MCMC sam- 
pling to explore this degeneracy, we find that the maximum likeli- 
hood spectral indices are a = -0. 1 1 ± 0.09 and a' = 2.32 ± 0.04, 
where (as previously in this work) a and a' represent the spec- 
tral indices of power-laws which dominate at low frequency and 
high frequency, respectively. A spectral index of a = —0. 1 1 is typ- 
ical of optically thin free-free emission. Melis et al. (2011) define 
IRAS 04368-1-2557 as a Class object and a free-free component 
would not be unusual for Class objects which predominantly 
have radio counterparts. We find an opacity index /3' ~ a' — 2 = 
0.32 ± 0.04 substantially shallower than that of Melis et al. (201 1), 
and consistent with the presence of large dust grains, or pebbles, in 
the disk. 



4 DISCUSSION 
4.1 Spectral Indices 

Using sub-millimetre data to derive both opacity indices and disk 
masses is common practice, however these values are uncertain 
due to both the approximately linear dependency on dust temper- 
ature, Td, and also the potentially significant fraction of the sub- 
millimetre emission that is optically thick. Perhaps more impor- 
tantly, at shorter wavelengths the measured emission cannot be ex- 
clusively attributed to the disk alone and may be contaminated by 
the protostellar envelope. Measuring the emission at longer wave- 
lengths reduces the effect of at least two of these issues significantly 
as the optically thick portion of the emission will be significantly 
reduced, if not negligible, by centimetre wavelengths, and addition- 
ally at frequencies below ~ 300 GHz the disk component will be 
expected to dominate the SED in the majority of cases. 

Consequently, although the collisional growth of dust grains in 
the disks around young stars has been inferred from sub-millimetre 
observations, it is preferable to measure its properties directly from 
longer wavelength observations. The theoretically determined val- 
ues of the opacity index for populations of these large grains (j8 ^ 
I; Miyake & Nakagawa 1993; Dullemond & Dominik 2005) are 
consistent with those those measured from the centimetre spectrum 
of the objects in this sample. In the majority of cases for this sam- 
ple, the contribution of additional radio emission mechanisms at 
1.8 cm has found to be sub-dominant with ^ 70 per cent of the flux 
density attributable to thermal dust emission. IRAS 04368-1-2557 is 
an exception to this with only 20 per cent of the 1.8 cm emission 
coming from thermal dust. This is notable as IRAS 04368-1-2557 is 
the only object in this sample which is Class 0, rather than Class I. 
This lack of free-free emission in the wider sample is of particular 
note as such emission in young stellar objects is often attributed 
to the presence of a molecular outflow, with shock ionization in 
the outflow proposed to be the source of the ionised plasma giving 
rise to the free-free emission. All the objects within this sample 
are known to possess molecular outflows (Gramajo et al. 2010), 
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IRAS 04368-1-2557: 




4''40"M0' 40"'00=' 39"50' 39'"40= 

Right Ascension (J2000) 




1 00 1 000 1 0^ 

Frequency [OHi] 



Figure 5. Top: AMI-LA 1.8 cm map of IRAS 04368-1-2557. Contours are 
shown in increments of a from 5 a, where a for is listed in Table 1 . The po- 
sition of the target object is indicated by a white cross ('-I-'), the FWHM of 
the AMI-LA primary beam by a circle, and the synthesized beam is shown 
as a filled ellipse in the bottom left comer of the map. This map is uncor- 
rected for primary beam attenuation. Bottom: Spectral energy distribution 
for IRAS 04368-H2557. Data points at V > 200 GHz are taken from Gra- 
majo et al. (2010), data points at lower frequencies are compiled from the 
literature (see text for details) and this work, see Table 2. A best fitting SED 
from the grid models of Robitaille et al. (2007) is shown as a solid line at 
V > 200 GHz, with the disk component of this model shown as a dashed 
line. Below v = 300 GHz this disk component is extrapolated as a power- 
law, see text for details, and shown as a dot-dash line. A power-law radio 
spectrum is also shown fitted to lower frequency radio data at v < 10 GHz 
with a = —0. 1 1, see text for details. 



however as has been observed in larger samples of Class & I 
objects the detection rate at cm-wavelengths for Class I sources 
(~ 20 per cent) is significantly lower than that of Class sources 
(~ 70 per cent; AMI Consortium: Scaife et al. 201 1). However we 
note that the contribution of both thermal and non-thermal emis- 
sion in young stellar objects appears to be variable (e.g. Shirley 
et al. 2007) as indicated by the discrepant flux densities at radio 
frequencies from observations at different epochs in the literature. 
Extrapolating power-law SEDs from the sub-mm to the AMI- 
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IRAS 04325+2402: 



IRAS 04361+2547: 
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Figure 4. Top: AMI-LA 1.8 cm maps of (a) IRAS 04325+2402b & (b) IRAS 04361+2547. Contours are shown in increments of a from 5 (T, where a for 
each object is listed in Table 1. The position of the target object is indicated by a white cross ('+'), the FWHM of the AMI-LA primary beam by a circle, and 
the synthesized beam is shown as a filled ellipse in the bottom left comer of each map. These maps are uncorrected for primary beam attenuation. Bottom: 
Spectral energy distributions for (a) IRAS 04325+2402b & (b) IRAS 04361+2547. Data points at V > 200 GHz are taken from Gramajo et al. (2010), data 
points at lower frequencies are compiled from the literature (see text for details) and this work, see Table 2, with data from the AMI-LA indicated in red in 
the online version. A best fitting SED from the grid models of Robitaille et al. (2007) is shown as a solid line at v > 200 GHz, with the disk component of this 
model shown as a dashed line. Below V = 300 GHz this disk component is extrapolated as a power-law, see text for details, and shown as a dot-dash line. In 
the case of IRAS 04361+2547 a second power-law spectrum is also shown fitted to lower frequency radio data at V < 10 GHz with a = — 3.5. The combined 
spectrum from the two power-law components is shown as a solid line at V < 300 GHz. 



LA waveband requires shallow spectral indices of a = 2.0 — 2.7, 
with an average value of a = 2.26 corresponding to /3 = 0.26 ± 
0.22. These indices are consistent with grain growth in these ob- 
jects providing a population of larger dust grains or pebbles. Such 
dust populations, with grains exceeding centimetre sizes (Draine 
2005) have been theoretically demonstrated to have opacity indices 
approaching zero (Draine 2005; Miyake & Nakagawa 1993; Dulle- 
mond & Dominik 2005) and the SEDs for the sample presented in 
this work support these predictions. 



4.2 Disk Masses 

Using the SED fits to the radio and sub-mm data we can separate 
the contributions from high frequency dust emission and low fre- 
quency radio emission at 1.8 cm. The contribution of the thermal 
dust emission at this wavelength is listed for each source in Table 3. 
Making the assumption that all dust emission at this wavelength 
will be optically thin we calculate disk masses using 

We use a dust temperature of = 20 K, the characteris- 
tic temperature found by Andrews Sl Williams (2005) at 1.3 mm. 
We note that increasing the value of to that of 60 K assumed 
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Table 3. Opacity indices and disk masses for the sample. Column [1] lists the target name; column [2] the spectral index measured from AMI-LA data alone; 
colvmm [3] the spectral index dominant at low frequencies from the joint power-law fit; colimm[4] the spectral index dominant at higher frequencies from the 
joint power-law fit; column [5] the opacity index corresponding to the high frequency power-law; column [6] the flux density predicted by the high frequency 
power-law at 1.8 cm; column [7] the disk mass calculated from the flux density in column [6] assuming the opacity index J3'; column [7] the disk mass 
calculated from the flux density in colunm [6] assuming an opacity index of J3 = 1. 



Name Kami a a' j3' Spred.ie P = P = 1 

(/iJy) (Mq) (Mq) 



IRAS04016-H2610 


0.65 


-1.1 


2.4 


0.4 


399 


0.013 


0.150 


DGTauB 


1.79 


-0.1 


2.0 


0.0 


442 


0.003 


0.166 


IRAS04248+2612 






2.7 


0.7 


< 125 


< 0.012 


< 0.047 


IRAS04302+2247 




-0.1 


2.3 


0.3 


181 


0.004 


0.068 


IRAS04325+24{)2 


1.49 




2.1 


0.1 


238 


0.002 


0.089 


IRAS0436 1+2547 


2.82 


-3.5 


2.0 


0.0 


449 


0.003 


0.168 


IRAS04368+2557 


1.17 


-0.1 


2.3 


0.3 


174 


0.004 


0.065 



by Terebey et al. (1992) at 2.7 mm would decrease our derived 
mass values by a factor of ~ 3. As well as temperature, the disk 
mass is a sensitive function of jS due to its linear dependence 
on the opacity, xv. Following Beckwith et al. (1990), we assume 
Kv = 0.1(v/1000GHz)'^ cm^g-i (a dust to gas ratio of 1:100). 
With J3 = 1 this leads to a value for the opacity at 1.8 cm (16 GHz) 
of x:i6GHz = 0.002 cm2g-i. 

The disk mass values derived from sub-mm data by Andrews 
& Williams (2005) assume a value of J3 = 1 when calculating the 
opacity. Disk mass estimates made directly from the AMI-LA flux 
densities under the same assumption are systematically higher than 
those derived from sub-mm data, see Table 3. Drawing a compari- 
son with the disk mass values from the best fitting SED models of 
Robitaille et al. (2007) plotted in Figs. 1, 2 & 5 for each source, 
there is a loose agreement between the sub-mm disk mass esti- 
mates, although the scatter is large. However, the cm- wave mass es- 
timates from the AMI-LA data are systematically higher than those 
derived from shorter wavelengths. One possible physical explana- 
tion for this difference is that the larger mass estimates arise as 
consequence of increased emission from larger dust grains at cm- 
wavelengths through the disk. This emission does not contribute to 
the sub-mm flux density, which predominantly traces smaller grains 
nearer to the disk surface due to optical depth effects in the disk 
core. 

However, the disk mass estimates from the AMI-LA data us- 
ing values of J3 = J3' derived from the SED fits, see Table 3, are 
systematically smaller than those found in Andrews & Williams 
(2005). This is a strong function of the lower values of the opacity 
index and therefore the correspondingly larger opacities. A clear 
example of this is IRAS04361-I-2547, where the difference in calcu- 
lated disk mass is over an order of magnitude. It is imclear whether 
the apparent difference between the sub-mm derived mass values 
and those from the AMI-LA data is real. The sub-mm values use a 
value of J3 = 1 which is significantly larger than that derived from 
SED fitting in AW05, where j3 is found to have an average value 
of approximately zero, and so the sub-mm masses could be consid- 
ered to be overestimated. Conversely the opacity index is expected 
to flatten towards lower frequencies and therefore the disk masses 
derived at cm- wavelengths should probably assume lower values of 

Similarly to Andrews & Williams (2005) wc find that in the 
disk masses derived from the cm-wave flux densities assuming 
J3 = 1, although larger than those inferred from sub-mm data, are 
still often significantly lower than those required by theory for 



giant planet formation (Pollack et al. 1996; Boss 1998). Three 
objects are exceptions to this: IRAS040 16-1-26 10, DGTau B and 
IRAS04361-H2547. For each of these objects the calculated disk 
mass is above the limit required for planet formation of 0.13 M© 
(Boss 1998). 

The gravitational instability mechanism of giant planet for- 
mation is a short lived process and most likely to occur during 
the early protostellar stages, at a stellar lifetime of ^ 10^ yr (Boss 
1998). This lifetime is typical of Class I objects, and given the 
frequency of such planets observed around more evolved stellar 
objects, we might consequently expect that objects such as those 
found in this sample would in the majority have observationally 
measured masses compatible with such a process. As suggested by 
Andrews & Williams (2005), it is possible that an underestimation 
of the disk mass from the disk SED may arise from the presence 
of a population of large grains which do not contribute to the emis- 
sion at sub-mm to cm-wavelengths, but which acts to increase the 
mass of the disk. Such a situation is more likely at sub-mm than 
cm-wavelengths, however it is difficult to state that disk mass es- 
timates from sub-mm flux densities are universally underestimated 
as there will also be a competing contribution at these wavelengths 
from the protostellar envelope. This contribution can only be sep- 
arated out using high resolution interferometric observations. Al- 
ternatively, the core accretion mechanism of giant planet formation 
(e.g. Pollack et al. 1996) requires significantly longer timescales of 
order 10^ yr, allowing young stellar objects in the relatively early 
stages of their evolution to continue accreting mass onto their pro- 
tostellar disks and consequently to increase their mass. 



5 CONCLUSIONS 

We have extended the SEDs of seven young stellar objects down 

to cm-wavelengths in order to constrain that part of the spectrum 
dominated by their circumstellar disks. The key results from this 
sample are summarised here: 

• We find that for the Class I objects in this sample the emis- 
sion at 1.8 cm is still doirunated by the tail of the dust greybody 
spectrum, rather than alternative radio emission mechanisms, with 
70 — 100 per cent of the cm- wave emission being attributable to 
thermal dust emission. The same is not true for IRAS 043684-2557, 
which is assumed to be Class and has only a 20 per cent contri- 
bution from thermal dust at 1.8 cm. 
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• We find spectral indices consistent witti a flattening of tiie 
opacity index towards longer wavelengths, with an average value 
^ = 0.26 ± 0.22, consistent with a significant population of large 
dust grains, and in accordance with theoretical models for the col- 
lisional growth of grains within circumstellar disks. 

• We derive disk masses directly from the thermal dust contribu- 
tion to the cm- wave flux densities. Under an assumption of J3 = 1, 
following Andrews & Williams (2005) we find that these mass es- 
timates are systematically higher than those determined from sub- 
mm data. Since the average value of the opacity index from SED 
fitting is similar at 1.8 cm and 850 /xm, we attribute this differ- 
ence in disk masses to an increased emission contribution from 
larger dust grains within the disk at centimetre wavelengths. Un- 
der these assumptions disk masses in excess of the lower limit re- 
quired for giant planet formation (Boss 1998) are recovered in al- 
most 50 per cent of cases. 
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